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The  investigation  reported  herein  was  authorized  by  the  Los 
Angeles  Port  Authority  under  a contract.  Agreement  No.  WES  76-5,  dated 
12  December  1975- 

This  study  was  conducted  during  the  period  December  1975  to  April 
1976  in  the  Hydraulics  Laboratory,  U.  S.  Army  Engineer  Waterways  Experi 
ment  Station  (WES),  under  the  direction  of  Mr.  H.  B.  Simmons,  Chief  of 
the  Hydraulics  Laboratory,  Dr.  R.  W.  Whalin,  Chief  of  the  Wave  Dynamics 
Division  (WDD),  and  Mr.  C.  E.  Chatham,  Chief  of  the  Harbor  Wave  Action 
Branch.  Mr.  J.  R.  Houston  conducted  the  study  and  prepared  this  report 

Messrs.  R.  R.  Bottin,  C.  Bryant,  and  A.  W.  Garcia  (WDD)  prepared 
the  finite  element  numerical  grids.  Mr.  D.  G.  Outlaw  (WDD)  developed 
the  computer  plotting  programs  and  aided  in  eliminating  errors  in  the 
finite  element  grids.  Mr.  R.  R.  Bottin  plotted  the  frequency  response 
curves.  Dr.  D.  C.  Raney  assisted  in  running  the  computer  program  for 
several  grids. 

Mr.  H.  L.  Butler  (WDD)  rewrote  the  matrix-solving  subroutine  of 
the  hybrid  finite  element  computer  program,  thus  reducing  significantly 
the  running  time  of  the  program. 

A significant  proportion  of  the  numerical  computations  were  per- 
formed on  a CDC-7600  computer  at  the  Los  Alamos  Scientific  Laboratory 
(LASL),  Los  Alamos,  New  Mexico.  Mr.  P.  Johnson  and  Mrs.  Margaret 
McCormick  of  LASL  assisted  significantly  in  coordinating  visits  to  Los 
Alamos  and  scheduling  computer  time  on  the  LASL  system. 

Drs.  H.  S.  Chen  and  C.  C.  Mei  of  the  Massachusetts  Institute  of 
Technology  provided  documentation  of  the  hybrid  finite  element  computer 
program  they  developed  and  materials  to  aid  in  its  utilization. 

Directors  of  WES  during  the  conduct  of  the  study  and  the  prepara- 
tion and  publication  of  this  report  were  COL  G.  H.  Hilt,  CE,  and 
COL  John  L.  Cannon,  CE.  Technical  Director  was  Mr.  F.  R.  Brown. 
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CONVERSION  FACTORS,  U.  S.  CUSTOMARY  TO  METRIC  (Si) 

UNITS  OF  MEASUREMENT 

U.  S.  customary  units  of  measurement  used  in  this  report  can  be  con- 
verted to  metric  (Si)  units  as  follows: 

Multiply By  To  Obtain 

feet  0.30U8  metres 

feet  per  second  0.30E8  metres  per  second 

feet  per  second  per  second  0.30U8  metres  per  second  per  second 

square  feet  per  second  0.0929030*+  square  metres  per  second 

degrees  (angle)  0.017^+5329  radians 
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LOS  ANGELES  HARBOR  NUMERICAL  ANALYSIS 
OF  HARBOR  OSCILLATIONS 


PART  I:  INTRODUCTION 


Harbor  Oscillation  Problem 

1.  A harbor  excited  by  long  waves  with  a period  close  to  the 
natural  oscillation  period  of  the  harbor  experiences  a resonant  reaction 
that  may  produce  wave  heights  in  the  harbor  far  greater  than  the  inci- 
dent wave  heights.  Generally,  the  vertical  movement  of  water  in  the 
harbor  is  still  only  a few  tenths  of  a foot  due  to  the  very  small 
heights  of  incident  long  waves.  However,  the  horizontal  movement  of 
water  may  cover  a distance  greater  than  10  ft*  as  a result  of  the  long 
wavelength  of  the  incident  waves.  Such  horizontal  water  movements  cause 
surge  and  sway  motion  of  moored  ships  that  may  hamper  unloading  opera- 
tions and  damage  mooring  facilities. 

2.  A number  of  harbors  on  the  west  coast  of  the  United  States 
have  histories  of  ship  motion  problems  resulting  from  harbor  resonant 
oscillations.  It  is  likely  that  the  long-period  waves  responsible  for 
exciting  harbor  oscillations  permeate  the  Pacific  Basin.  Lack  of  modern 
harbor  development  may  prevent  the  occurrence  of  any  problem  for  most 
coastal  areas  of  the  Pacific.  Many  places  that  have  a problem  may 
merely  tolerate  it  through  ignorance  of  possible  remedial  measures  with 
the  result  that  the  problem  goes  unreported. 

3.  The  generating  mechanism  of  the  very  small  amplitude  long 
waves  that  excite  harbors  is  still  unknown.  However,  Munk's  surf-beat 
theory'1'  suggested  an  attractive  explanation  of  the  origin  of  these 
waves.  Surf  beats  result  from  the  superposition  of  two  trains  of  waves 


* A table  of  factors  for  converting  U.  S.  customary  units  of  measure- 
ment to  metric  (SI)  units  is  presented  on  page  3. 
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of  slightly  different  period.  The  beats  are  an  amplitude  modulation 
of  a shorter  period  wave  train.  Swell,  for  example,  almost  always  ex- 
hibits an  irregular  beating  up  to  several  minutes  in  period. 

4.  Carr  argued  that  local  surf  beat  could  not  cause  the  long- 
period  waves  known  to  excite  harbors.  Although  the  surf-beat  process 
produces  a long -period  vertical  water  motion,  the  really  important  fac- 

• tor  of  long-period  horizontal  water  motion  is  lacking  in  surf  beats. 

Severe  harbor  surging  also  is  known  to  occur  during  otherwise  calm  seas 
with  identical  long-period  wave  activity  reported  along  hundreds  of 
miles  of  a coastline.  Carr  maintained  that  such  activity  is  not  compat- 
ible with  long-wave  generation  by  local  surf  beat. 

5.  True  long-period  wave  trains  may  be  generated  from  surf  beats 
through  nonlinear  processes  in  the  surf  zone.  Munk"*'  estimated  that 
beaches  reflected  1 percent  of  the  energy  incident  upon  them  in  the  form 
of  these  long-period  waves.  It  is  perhaps  these  waves  reflected  from 
distant  shores  of  the  Pacific  that  cause  oscillation  problems  in  harbors 
on  the  west  coast  of  the  United  States.  Long-period  wave  activity  would 
then  increase  after  distant  storms  of  the  south  or  west  Pacific  Ocean 
created  large  swell  conditions  and  the  resulting  surf  beat  generated 
long-period  wave  trains  in  the  surf  zone  which  reflected  off  shorelines. 

3 

Wilson  has  noted  that  harbor  oscillation  in  Los  Angeles  and  Long  Beach 
Harbors  has  a seasonal  fluctuation  of  frequency  and  intensity  which  is 
! greatest  during  the  summer  months  and  least  during  the  winter.  He  noted 

$ that  such  a pattern  is  almost  identical  with  that  found  for  Table  Bay 

Harbour,  Cape  Town,  South  Africa,  and  suggests  that  the  storms  in  the 
Southern  Hemisphere  are  responsible  for  the  long-period  waves. 

6.  Further  evidence  that  long-period  wave  trains  exist  from  time 

to  time  as  a general  ocean  condition  and  are  not  generated  by  local 

L 

surf  beats  comes  from  a harbor  study  for  San  Nicholas  Bay,  Peru.  It 
was  concluded  from  horizontal  water  velocity  measurements  that  long- 
period  waves  affecting  San  Nicholas  Harbor  were  progressive  waves  un- 
associated with  the  concurrent  swell  waves. 
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Purpose  of  Study 

7.  The  purpose  of  this  study  was  to  investigate  harbor  oscilla- 
tions excited  by  long  waves  with  periods  from  1 to  10  min  using  a finite 
element  numerical  model.  Existing  conditions  of  Los  Angeles  and  Long 
Beach  Harbors  and  conditions  of  a proposed  harbor  expansion  were  con- 

i sidered  to  ensure  that  the  expansion  would  not  cause  undesirable  oscil- 

lations in  existing  basins  or  characteristic  oscillations  of  its  own 
that  were  undesirable. 

8.  The  amplification  peaks  predicted  by  the  numerical  model  may 
be  much  larger  than  the  peaks  which  actually  occur  in  nature,  since  the 
model  neglects  all  dissipative  processes  except  radiation  of  energy  by 
a harbor.  However,  the  model  adequately  predicts  the  relative  severity 
of  various  modes  of  oscillation,  and.  the  results  can  be  used  to  compare 
oscillation  characteristics  of  the  existing  Los  Angeles  and  Long  Beach 
Harbors  complex  and  proposed  harbor  expansions.  Information  about  pe- 
riods of  resonant  peaks  and  the  spatial  patterns  of  wave  heights  and 
currents  when  areas  of  the  harbor  are  excited  also  will  be  used  to  de- 
termine both  the  incident  wave  period  ranges  which  need  to  be  tested 
and  the  optimum  gage  placement  during  future  wave  tests  in  the  Los 
Angeles  and  Long  Beach  Harbors  hydraulic  model. 
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PART  II:  NUMERICAL  MODEL 


Response  of  the  Los  Angeles  and  Long  Beach  Harbors  complex  to 
excitation  was  determined  by  using  a hybrid  finite  element 
model  developed  recently  by  Chen  and  Mei  at  the  Massachusetts 
of  Technology . ^ The  model  solves  the  following  generalized 
equation: 

2 

V • [h(x,y)V<|)(x,y ) ] + <t>(x,y)  = 0 (l) 

§ 

where  <J>(x,y)*  is  the  velocity  potential  defined  by  U(x,y)  = -V<}>(x,y), 
with  U(x,y)  being  a two-dimensional  velocity  vector  and  w an  angu- 
lar frequency,  h(x,y)  is  the  water  depth,  and  g is  the  acceleration 
due  to  gravity.  Equation  1 governs  small  amplitude  undamped  oscilla- 
tions of  water  in  a basin  of  arbitrary  shape  and  variable  depth  forced 
by  periodic  l'jng  waves.  It  has  been  further  assumed  that  the  flow  is 
irrotatxonal . 

1U.  The  Helmholtz  equation: 

2 

V2<j>(x,y)  + ^ 4>(x,y)  = 0 (2) 

is  the  governing  equation  for  a constant-depth  ocean  region  outside  the 
basin. 

11.  For  a harbor  in  a semi-infinite  ocean  with  a straight  coast- 
line there  is  an  incident,  reflected,  and  scattered  wave.  The  scattered 
wave  has  a velocity  potential  <j>  given  by 

00 

<j>  = T a H (kr)  cos  n0  (3) 

s **  n n 
n=0 

where  a are  unknown  coefficients  and  H (kr)  are  Hankel  functions  of 
n n 

the  first  kind  of  order  n . 


9. 

long  wave 
numerical 
Institute 
Helmholtz 


* For  convenience,  symbols  and  unusual  abbreviations  are  listed  and 

defined  in  the  Notation:  Appendix  A. 


12.  satisfies  the  radiation  condition  that  the  scattered 

wave  must  behave  as  an  outgoing  wave  at  infinity.  This  condition  is 
known  as  the  Sommerfeld  radiation  condition  and  may  be  expressed 
mathematically  as  follows: 


r-w°  \ \3r 


0 


(U) 


13.  Chen  and  Mei  used  a calculus  of  variations  approach  and 
obtained  a Euler-Lagrange  formulation  of  the  boundary  value  problem. 

The  following  functional  with  the  property  that  it  is  stationary  with 
respect  to  arbitrary  first  variations  of  <J)(x,y)  was  constructed  by 
Chen  and  Mei: 

2 

F(4>)  = ff  1/2  {h(V<}>)2  - — <p2}  dA 

6 

3 ( 4>r  — 4>x  ^ 3 ( — 4>T ) } 

+ 1/2  §{h(<t  - O — da  - § {hi})  — da 

R I 3n  a 3n 

a a 

3 ( 4>R  - 4>T  ) 1 

- § {hi})  -—■}  da  + §{hi}>T  — r da  (5) 

a 3n  I 3n 

a a 


where 

A = the  region  inside  the  harbor 

§ = line  integral 

4>  = far  field  velocity  potential 

n 

<J> = velocity  potential  of  the  incident  wave 

n&  = unit  normal  vector  outward  from  region  A 
a = boundary  of  region  A 

<}i  = total  velocity  potential  evaluated  on  boundary  a 

3. 

lh.  Proof  was  given  by  Chen  and  Mei  that  the  stationarity  of 
this  functional  is  equivalent  to  the  original  boundary  value  problem. 

15.  The  integral  equation  obtained  from  extremizing  the  func- 
tional is  solved  by  utilizing  the  finite  element  method.  This  method 
is  a technique  of  numerical  approximation  that  involves  dividing  a 
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domain  into  a number  of  nonoverlapping  subdomains  which  are  called 
elements . 

16.  The  solution  of  the  problem  is  approximated  within  each 
element  by  suitable  interpolation  functions  in  terms  of  a finite  number 
of  unknown  parameters.  These  unknown  parameters  are  the  values  of  the 
field  variable  <p(x,y)  at  a finite  number  of  points  which  are  called 
nodes.  The  relations  for  individual  elements  are  combined  into  a 
system  of  equations  for  all  unknown  parameters. 

IT.  In  the  region  outside  the  basin,  the  velocity  potentials 
are  solved  analytically  in  terms  of  unknown  coefficients.  The  region 
is  considered  a single  element  with  an  "interpolation  function"  given 
by  Equation  3.  The  infinite  series  is  terminated  at  some  finite  value 
such  that  the  addition  of  further .terms  does  not  significantly  influence 
the  calculated  values  of  <}>(x,y)  . The  resulting  equation  is  combined 
with  the  system  of  equations  for  unknown  parameters  at  nodal  points 
within  the  basin  and  this  complete  system  is  solved  using  Gaussian 
elimination  matrix  methods. 

18.  n(x,y)  is  related  to  <f>(x,y)  through  the  linearized  dynamic 
free  surface  boundary  condition 


n(x,y ) 


1.  94>(x,y ) 
g 3t 


(6) 
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19-  The  horizontal  velocity  components  have  the  following  form: 


u(x,y)  - - a Mixii  . v(XjJ)  . as 

w 3x  w 3y 


(T) 


20.  The  hybrid  finite  element  method  (so  named  by  Chen  and  Mei 
because  the  method  involves  the  combination  of  analytical  and  finite 
element  numerical  solutions)  is  a steady-state  solution  of  the  boundary 
value  problem.  The  response  of  a harbor  to  an  arbitrary  forcing  func- 
tion can  be  easily  determined  within  the  framework  of  a linearized 
theory.  For  example,  an  arbitrary  incident  wave  amplitude  at  the 
harbor  mouth  in  the  absence  of  the  harbor  can  be  Fourier  decomposed  as 
follows: 


r 

t 

i 
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dw 


b (t) 

O 


b(w) 


-iwt 

e 


If  n(x,y, w)  is  the  response  amplitude  at  any  point  (x,y)  inside  the 
harbor  due  to  an  incident  plane  wave  of  unit  amplitude  and  frequency 
w , then  the  response  of  the  harbor  to  the  arbitrary  incident  wave 
amplitude  bQ(t)  is  given  by 


£(x,y,t) 


b(w)n(x,y ,w)e  iwt 


dw 


(9) 


where  the  operation  R [ ] takes  the  real  part  of  the  bracketed 

quantity. 

21.  Therefore,  as  soon  as  n(x,y,w)  is  known  for  all  w , the 
harbor  response  to  an  arbitrary  incident  wave  can  be  calculated. 

22.  Chen  and  Mei^  verified  their  finite  element  numerical  model 
by  comparing  the  model's  calculations  with  analytical  and  experimental 
results  for  simple  wave  problems.  Figure  1 is  a comparison  of  analyt- 
ical and  numerical  results  of  wave  diffraction  by  a vertical  circular 
cylinder.  Similar  excellent  agreement  was  obtained  by  comparing  numeri- 
cal calculations  with  an  analytical  solution  for  a semicircular  harbor 
with  two  openings.  Experimental  data  for  the  oscillation  of  a rectangu- 
lar harbor  also  were  compared  with  the  finite  element  numerical  calcu- 
lations. Calculated  periods  of  the  resonant  peaks  agreed  with  experi- 
mental data;  however,  calculated  amplitudes  of  the  resonant  peaks  were 
larger  since  the  numerical  model  did  not  consider  frictional  effects. 

23.  Several  modifications  were  made  to  the  numerical  computer 
program  of  Reference  5 before  calculations  were  made  for  Los  Angeles 
and  Long  Beach  Harbors.  Slight  modifications  were  made  to  allow  the 
model  to  handle  harbors  with  variable  depths.  Chen  and  Mei^  applied 
their  numerical  model  only  to  constant-depth  problems.  Hie  subroutine 
which  solves  the  large  system  of  algebraic  equations  arising  from 
finite  elements  was  modified  to  be  more  efficient.  This  modification 
resulted  in  a decrease  in  the  computational  time  of  the  numerical  model 
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Figure  1.  Comparison  of  numerical  and 
analytical  results  for  wave  scattering 
off  a vertical  circular  cylinder^ 
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by  a factor  of  between  two  and  three. 

2k.  Velocity  calculations  also  were  added  to  the  numerical 
model. ^ Velocities  are  calculated  using  the  expressions  in  Equation  7. 
The  slope  of  n may  be  obtained  at  any  point  (x,y)  within  an  element 
through  the  matrix  equation 


in 

3x 

in 

3y 


G{  n }e 


where  the  matrix  vector 


(nr 


is  defined  as 


(n) 


e 


(10) 


(11) 


11 


f 


and  the  coordinates  of  the  three  nodes  (i,j,k)  forming  the  elements  are 
defined  as  (x^y^),  (x^y^),  and  (x^yj^). 

25.  G is  the  element  slope  matrix  defined  as 


’3N. 

1 

3N. 

3N 

k 

3x 

9x 

3x 

3N. 

3N. 

3N 

1 

1 

k 

_3y 

3y 

3y 

and  matrix  N is  known  as  the  "interpolation"  function  and  is  defined 
as 


N = {N . ,N . ,N  } 
1 J k 


where 


Ni  = {(xkyj  ~ Xjyk}  + (yk  ' yj)x  + Uj  " Xk)y)/2A 

N . ,N  = similar  expressions  obtained  by  the  cyclical  permutation 
J of  i , j , and  k 

0 

A = area  of  the  element  A 


26.  The  {n}  matrix  is  known  from  the  finite  element  computa- 
tions. Together  with  the  interpolation  function  N , the  horizontal 
velocity  at  any  point  (x,y)  within  an  element  can  be  determined.  Ve- 
locities presented  in  this  report  were  calculated  at  element  centers. 

27.  The  finite  element  numerical  model  has  several  important 

advantages  compared  with  other  numerical  harbor  oscillation  models  such 

7 

as  that  of  Lee  and  Raichlen.  First,  the  model  allows  a variable  depth. 
Furthermore,  the  computational  time  requirement  for  the  model  is  sub- 
stantially less  than  computational  times  required  by  other  numerical 

g 

models.  It  is  the  only  present  numerical  model  for  harbor  oscillations 
which  can  feasibly  calculate  resonance  effects  in  a harbor  as  large  as 
the  Los  Angeles  and  Long  Beach  Harbors  complex. 
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PART  III:  APPLICATION  OF  NUMERICAL  MODEL 

Boundary  Conditions 

28.  The  boundary  condition  used  in  the  numerical  model  along  the 
shoreline  of  the  Los  Angeles  and  Long  Beach  Harbors  complex  is  that  the 
normal  component  of  the  velocity  be  equal  to  zero.  The  same  boundary 
condition  was  used  for  the  detached  breakwater  surrounding  the  proposed 
oil  tanker  terminal  at  Pier  J in  Long  Beach  Harbor.  The  breakwater  core 
extends  above  still  water  and  is  impermeable  to  incident  wave  energy  and 
is  therefore  represented  as  a solid  barrier  in  this  study.  The  eastern 
entrance  to  the  tanker  terminal  basin  is  600  ft  wide  and  the  basin  is 
dredged  to  a depth  of  62  ft  mllw. 

29.  The  San  Pedro  breakwater  (Plate  l)  does  not  have  a solid 
core.  It  is  composed  of  class  B rock  having  an  average  diameter  of  ap- 
proximately 1.8  ft  and  has  a crown  at  mllw  of  fitted  concrete  blocks. 

9 

Keulegan  gives  the  following  formula  for  calculating  wave  transmission 
through  rock  structures: 


(iM 

(15) 

and 

= incident  wave  amplitude 
a^  = transmitted  wave  amplitude 
h = water  depth 
L = width  of  breakwater 
X = wavelength  of  incident  waves 
T = wave  period 
Using  these  parameters 
h = 50  ft 


a. 

a^  1 + ^ V h / X 


(r) 


where 


U/3 


X / T 

y^.H^fgh- 
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r 


L = 125  ft 
T = 180  sec 
X = 7222  ft 
d = 1.8  ft 


ai  = 0.05  ft 


yields 


— = 0.87 

3, . 


30.  This  figure  is  a lower  bound  for  the  transmission  through  the 
San  Pedro  breakwater  for  a 3-min  wave  since  a^  = 0.05  ft  is  probably 
the  upper  bound  wave  amplitude  for  a 3-min  wave  at  the  San  Pedro  break- 
water. This  upper  bound  can  be  deduced  from  prototype  measurements 
taken  in  Los  Angeles  and  Long  Beach  Harbors  by  the  WES."^  Long-period 
wave  activity  cannot  be  discerned  in  wave  gage  measurements  at  the  San 
Pedro  breakwater  because  of  the  large  wave  heights  in  the  swell  regime. 
However,  a gage  off  Reservation  Point  in  the  Los  Angeles  Harbor  was 
sufficiently  protected  from  swell  that  visual  inspections  could  be  made 
of  long-period  wave  activity.  A total  of  107  visual  inspections  of 
4-hr  analog  charts  taken  during  ship  motion  occurrences  from  a year's 
prototype  recording  reveal  only  a single  instance  of  a wave  height  as 
great  as  0.1  ft  for  waves  with  periods  of  1 to  10  min.  The  single  re- 
cording was  of  a 6-min  wave  and  probably  results  from  the  large  6-min 
resonant  oscillation  that  occurs  in  the  East  Channel  of  Los  Angeles 
Harbor . 

31.  A 3-min  wave  was  considered  in  the  above  calculations  since 
longer  period  waves  of  the  same  amplitude  would  more  readily  be  trans- 
mitted through  the  breakwater  and  the  finite  element  grids  for  1-  to 
3-min  waves  (paragraphs  3*+-36)  do  not  extend  as  far  as  this  breakwater. 

32.  Since  the  San  Pedro  breakwater  is  very  pervious  to  3-  to 
10-min  incident  waves,  it  is  not  represented  in  the  finite  element 
grids.  Although  the  breakwater  looks  impervious  from  the  water  as  a re- 
sult of  the  fitted  concrete  blocks  extending  up  from  the  waterline,  Wil- 

3 „ 

son  noted  that  this  breakwater  has  long  been  considered  an  ineffective 
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filter  against  long  waves  because  of  its  1*3  percent  voids;  this  is  in- 
deed the  case. 

33-  The  Middle  breakwater  (Plate  l)  of  Los  Angeles  and  Long 
Beach  Harbors  has  a solid  core  to  -26  ft  mllw  which  is  covered  with  a 
layer  of  class  B rock  and  a layer  of  class  A rock  (approximately  3 ft 
in  diameter).  The  rock  covering  of  the  solid  core  is  at  least  as  per- 
meable to  long  waves  as  the  San  Pedro  breakwater,  since  it  consists  of 
rock  the  same  size  or  larger.  Therefore,  the  Middle  breakwater  is  rep- 
resented by  a depth  change  with  the  core  impermeable,  and  small  finite 
element  triangles  are  placed  along  the  breakwater  in  the  numerical  grid 
to  properly  represent  it. 


Numerical  Grids 


3^+.  The  numerical  grids  used  in  the  numerical  model  computations 

are  shown  in  Plates  3-12.  The  elements  are  all  triangular  and  often 

equilateral.  The  lengths  of  the  triangle  sides  were  in  general  selected 

so  that  they  were  one  eighth  of  the  local  wavelength  or  less.  Chen  and 

Mei^  compared  the  finite  element  solution  of  diffraction  by  a vertical 

cylinder  with  an  exact  analytical  solution  and  found  a global  error  of 

0.7  percent  when  they  used  an  element  length  of  one  tenth  the  incident 

1 E 

wavelength.  They  defined  a global  error  as  — (local  error),  where 
the  summation  is  performed  over  n finite  element  nodes.  The  local 
error  is  defined  as 


n - nt 

Local  error  = (l6) 

\ 

where  p is  the  magnitude  calculated  by  the  finite  element  method  at 
a nodal  point  and  p^  is  the  theoretical  value.  In  the  present  study, 
it  was  judged  that  element  lengths  equal  to  one  eighth  of  the  local 
wavelength  or  less  would  provide  sufficient  accuracy.  Element  lengths 
vary  in  the  grids  in  relation  to  local  depths.  Element  lengths  near 
harbor  entrances  are  made  much  smaller  than  may  be  strictly  necessary 
in  order  to  ensure  proper  propagation  of  waves  into  inner  basins. 
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35.  Five  grids  were  used  for  each  case  tested  and  are  tabulated 
below.  It  is  not  feasible  to  use  a single  grid  to  calculate  harbor  re- 
sponse for  incident  waves  with  periods  from  1 to  10  min.  For  example, 

1-  to  3-min  grids  could  not  be  constructed  to  cover  the  entire  harbor 
complex  because  of  computer  memory  and  running  time  limitations.  Sepa- 
rate grids  were  used  for  3-  to  U-,  U-  to  6-,  and  6-  to  10-min  calcula- 
tions to  reduce  the  cost  of  operating  the  model.  The  computational  time 
of  the  model  is  approximately  proportional  to  the  number  of  nodes  in  the 
grid  times  the  bandwidth  squared.  The  bandwidth  is  the  greater  of 
either  the  maximum  difference  for  all  elements  between  two  adjacent  node 
numbers  in  an  element  of  the  grid  or  the  sum  of  the  number  of  nodes  on 
the  grid  semicircle  plus  the  number  of  terms  retained  in  truncating  the 
infinite  series  of  Equation  3.  In  general,  the  bandwidth  increases  ap- 
proximately linearly  with  a decrease  in  the  lengths  of  the  triangular 
element  sides.  The  number  of  nodes  increases  approximately  as  the  in- 
verse square  of  the  triangle  side  lengths.  Thus  a grid  for  3-min  waves 
(eight  nodes  per  wavelength)  requires  approximately  l6  times  more  com- 
putational time  than  one  for  6-min  waves  (eight  nodes  per  wavelength). 


Existing  Conditions  Revised  Conditions 


Grid 

1 

1 

to 

2 min 

Grid 

1A 

1 

to 

2 min 

Grid 

2 

2 

to 

3 min 

Grid 

2A 

2 

to 

3 min 

Grid 

3 

3 

to 

U min 

Grid 

3A 

3 

to 

It  min 

Grid 

k 

k 

to 

6 min 

Grid 

UA 

U 

to 

6 min 

Grid 

5 

6 

to 

10  min 

Grid 

5A 

6 

to 

10  min 

36.  Many  of  the  numerical  grids  for  existing  conditions  were  de- 
signed so  that  they  could  be  more  readily  modified  to  construct  grids 
for  the  harbor  expansion  conditions.  However,  the  design  of  the  harbor 
was  modified  several  times  during  the  course  of  this  study.  Therefore, 
some  element  patterns  in  the  grids  may  have  no  particular  purpose. 

Grids  also  were  constructed  to  minimize  the  bandwidth,  and  this  require- 
ment occasionally  governed  triangle  shapes  for  areas  of  the  grids. 
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PART  IV : RESULTS 


Comparison  of  Prototype  Data  and  Numerical 
Calculations  for  Existing  Conditions 

37-  Prototype  wave  gage  recordings  were  collected  at  Los  Angeles 
and  Long  Beach  Harbors  from  26  May  1971  to  29  June  1972  by  the  U.  S. 

Army  Engineer  Waterways  Experiment  Station.10  Energy  spectra  were  de- 
termined for  10  time  periods  of  the  year  during  which  ship  motion  prob- 
lems were  reported.  Spectral  analyses  were  performed  via  Fast  Fourier 
Transform  (FFT)  techniques.  Variability  of  spectral  estimates  was  de- 
creased by  smoothing  over  10  adjacent  frequency  bands.  Ten  bands  were 
chosen  so  that  there  would  still  be  sufficient  frequency  resolution  to 
define  the  long-period  energy  for  wave  periods  as  great  as  10  min.  The 
data  record  length  analyzed  for  each  spectral  analysis  was  21,576  sec 
(6.83  hr)  and  the  lag  time  between  successive  record  start  times  was 
chosen  to  be  3 hr.  Further  details  of  the  analyses  can  be  found  in 
Reference  10. 

38.  Tables  1 and  2 give  estimates  of  wave  periods  and  spectral 
energy  densities  of  peak  energies  for  prototype  wave  gages  1 and  5; 

Plate  1 shows  the  locations  of  these  gages.  The  starting  times  in  Ta- 
bles 1 and  2 give  the  year,  Julian  day,  and  hour  on  a 2l-hr  clock.  Each 
of  the  10  time  periods  lasts  approximately  2 days. 

39-  Tables  1 and  2 give  the  wave  period  with  the  greatest  spec- 
tral energy  density  within  a range  of  periods.  The  range  of  periods  is 
such  that  the  spectral  lines  within  the  range  cannot  be  resolved  by  the 
spectral  analysis;  for  example,  in  the  range  of  180-360  sec  the  adja- 
cent frequency  lines  corresponding  to  periods  of  159-36  and  387-02  sec 
cannot  be  resolved  by  the  spectral  analysis  and  should  be  considered  the 
same . 

10.  Plate  13  is  a time-dependent  spectral  plot  of  gage  1 data  for 
the  time  period  of  1500  November  27  through  0300  November  29,  1971-  The 
variability  of  energy  levels  in  this  plate  illustrates  the  nonstation- 
arity  of  the  variance  in  energy  bands.  Long-period  waves  are  in  general 
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transient  (time-dependent)  phenomena,  and  some  caution  must  be  observed 
in  the  application  of  spectral  analysis  techniques  to  nonstationary 
processes.  However,  Tables  1 and  2 do  give  an  indication  of  character- 
istic energy  bands  for  each  gage  location. 

41.  It  should  be  noted  that  the  finite  element  model  calculations 
provide  only  the  response  characteristics  of  Los  Angeles  and  Long  Beach 
Harbors.  This  response  must  be  multiplied  by  the  input  spectrum  of 
waves  for  a comparison  with  prototype  data.  However,  the  input  wave 
spectrum  is  not  known  since  long  waves  arriving  from  the  deep  ocean  had 
heights  too  small  to  be  measured  by  gages  near  the  outer  breakwaters 
during  the  prototype  collection  (see  paragraph  30).  Still,  a corre- 
spondence exists  between  resonant  peaks  calculated  by  the  numerical  model 
and  those  measured  in  the  prototype.  Plates  14-18  and  24-28  are  plots 

of  the  wave-height  amplification  factor  calculated  by  the  finite  ele- 
ment numerical  model  for  prototype  gages  1 and  5,  respectively.  The 
frequency  response  for  the  existing  conditions  of  Los  Angeles  and  Long 
Beach  Harbors  shown  in  Plate  1 is  compared  with  the  frequency  response 
for  the  revised  conditions  shown  in  Plate  2.  The  wave-height  amplifi- 
cation factor  is  defined  at  each  point  as  the  wave  height  at  the  point 
divided  by  twice  the  incident  wave  height.  This  definition  of  amplifi- 
cation factor  is  traditional  and  is  a result  of  the  fact  that  the  stand- 
ing wave  height  for  a straight  coast  with  no  harbor  would  be  twice  the 
incident  wave  height  due  to  the  superposition  of  the  incident  and  re- 
flected waves. 

42.  Table  1 shows  a large  resonant  peak  at  387  sec  for  gage  1. 
Plate  18  shows  large  wave-height  amplification  for  the  period  range  of 
360-480  sec.  As  was  mentioned  earlier,  the  adjacent  frequency  lines 
corresponding  to  periods  of  459*36  and  387*02  sec  in  the  range  360- 

480  sec  cannot  be  resolved  by  spectral  analysis  and  should  be  considered 
the  same.  Appendix  C of  Reference  1 presents  a spectral  analysis  of 
data  taken  over  a 2-day  period  for  gage  1 which  has  not  been  smoothed 
over  10  adjacent  frequency  bands  to  reduce  variability.  Resonant  peaks 
occur  at  482,  4l0,  372,  and  346  sec  in  this  analysis.  The  period  loca- 
tion of  these  peaks  is  in  good  agreement  with  Plates  17  and  18. 


Plate  16  shows  a moderately  large  resonant  peak  at  233  sec  and  spectral 
plots  in  Reference  10  show  a small  but  well-defined  peak  at  263  sec. 
Table  1 shows  a resonant  peak  from  115  to  ll2  sec  and  Plate  15  shows  a 
similar  peak  predicted  by  the  numerical  model.  Table  1 also  shows  a 
small  peak  from  66  to  78  sec.  Plate  ll  shows  a similar  peak,  but  it  is 
no  -larger  than  peaks  from  85  to  95  sec  and  105  to  115  sec  also  shown  in 
Plate  lit.  The  water  depths  in  the  dredge  spoil  area  just  outside  the 
East  Channel  which  were  used  on  the  numerical  model  were  taken  from 
National  Ocean  Survey  (formerly  U.  S.  Coast  and  Geodetic  Survey)  maps 
published  in  1973.  It  is  possible  that  depths  in  this  area  have  changed 
enough  since  the  prototype  measurement  to  influence  some  of  the  oscil- 
lation characteristics  of  the  East  Channel.  The  resonant  peaks  dis- 
played in  Table  1 are  generally  in  fair  agreement  with  the  numerical 
results  shown  in  Plates  ll-l8. 

1+3.  A comparison  of  predictions  of  the  numerical  model  and  spec- 
tral analysis  of  the  prototype  data  for  five  gages,  prototype  gages  5» 

6,  7,  8,  and  11,  was  made  in  Reference  11.  These  gages  were  located  in 
the  Long  Beach  Harbor  complex.  Excellent  agreement  was  found  between 
the  numerical  results  and  the  prototype  data  for  all  five  gages.  Ta- 
ble 2 shows  estimates  of  wave  periods  and  spectral  energy  density  of 
peak  energies  for  gage  5 and  Plates  21+-28  show  numerical  model  calcu- 
lations of  the  wave-height  amplification  factor  at  this  gage.  The  solid 
line  in  these  plots  shows  calculations  for  existing  conditions. 

Table  2 indicates  that  there  is  no  resonant  peak  for  gage  5 
data  between  8 and  10  min.  Numerical  calculations  presented  in  Plate  28 
also  show  no  resonant  peak  for  this  period  range.  Table  2 shows  large 
peaks  at  387  and  237  sec.  Plates  26-28  show  large  peaks  at  396  and 
222  sec  and  no  peaks  from  2l0  to  360  sec.  Neither  the  numerical  compu- 
tations in  Plate  25  nor  the  prototype  data  of  Table  2 indicate  peaks 
from  120  to  210  sec.  The  prototype  data  have  peaks  at  100,  93,  and 
61  sec.  The  numerical  model  calculates  peaks  at  92,  89,  and  63  sec 
(Plate  2l).  Excellent  agreement  on  the  resonant  peaks  is  apparent  be- 
tween prototype  data  and  the  numerical  model  calculations. 

15.  A detailed  analysis  of  the  prototype  data  for  gages  2 and  11 
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was  not  performed  in  Reference  10  because  no  significant  oscillations 
from  1 to  10  min  were  apparent.  Time-dependent  spectral  plots  for 
gages  2 and  11  in  Reference  10  indeed  show  no  significant  long-period 
oscillations.  Plates  18-23  show  numerical  model  calculations  of  the 
wave-height  amplification  factor  at  gage  2.  No  significant  resonant 
peaks  occur  from  1 to  10  min.  Similarly,  the  numerical  model  does  not 
predict  large  resonant  peaks  at  the  gage  11  location.  Plots  for  gage  11 
have  been  omitted. 
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Comparison  of  Harbor  Response  for  Existing 
and  Revised  Conditions 


hG.  Plates  lb-23  show  a comparison  of  plots  of  the  wave-height 
amplification  factor  as  a function  of  wave  period  at  gages  1 and  2 for 
existing  and  revised  harbor  conditions.  Plates  1 and  2 show  the  exist- 
ing harbor  and  the  proposed  modifications.  The  modifications  include 
landfills  and  dredge  cuts.  Wave  periods  from  1 to  10  min  are  considered. 

hj . Plates  29-33  are  frequency  response  plots  of  wave-height 
amplification  for  the  proposed  Liquified  Natural  Gas  (LNG)  slip  (sta  l) 
of  Los  Angeles  Harbor.  The  largest  resonant  peaks  occur  at  110,  2l6, 

238,  and  258  sec.  The  258-sec  peak  is  fairly  large  and  the  mode  of 
oscillation  rather  complicated.  Plate  187  shows  contours  of  the  wave- 
height  amplification  factor  for  this  resonant  peak.  The  oscillation  of 
the  LNG  slip  is  in  phase  with  an  oscillation  of  the  end  basin  of  Fish 
Harbor  and  180  degrees  out  of  phase  with  the  proposed  general  cargL 
terminal  to  be  created  by  the  Terminal  Island  landfill  and  the  center 
basin  of  Fish  Harbor.  These  harbor  areas  are  undergoing  coupled  oscil- 
lations with  water  flowing  back  and  forth  between  different  areas. 

UQ.  Plates  3b-b2  are  frequency  response  plots  of  wave-height  am- 
plification for  sites  around  the  Los  Angeles  Harbor  complex  at  which 
prototype  gages  were  not  located  but  which  may  experience  resonant 
problems  (LA-1  through  LA-5).  Plate  2 shows  locations  of  these  gages. 
Plots  of  the  frequency  response  for  gage  LA-2  at  the  northern  tip  of 
the  Watchorn  Basin  are  not  included  in  this  report  because  areas  of 
this  basin  were  too  shallow  to  allow  adequate  resolution  of  the  waves. 
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This  basin  is  a small-boat  harbor  and  waves  with  periods  from  1 to 
10  min  do  not  cause  ship  motion  problems  for  such  small  craft.  Fre- 
quency response  plots  of  wave-height  amplification  for  gages  LA-3  and 
LA- 4 are  not  shown  since  they  are  similar  to  the  plots  for  prototype 
gage  1 but  have  smaller  peaks.  Current  velocity  plots  for  these  two 
gages  are  presented  later. 

1+9.  Resonant  peaks  of  the  wave-height  amplification  factor  for 
gages  LA-1  and  LA-5  occur  at  about  the  same  wave  periods  for  existing 
and  modified  conditions.  Sometimes  the  resonant  peaks  are  larger  for 
one  of  these  conditions  and  sometimes  for  the  other.  Many  of  these  dif- 
ferences are  probably  attributable  to  difficulties  in  locating  the  peaks 
and  not  differences  in  response.  Response  calculations  were  initially 
made  every  2,  3,  and  5 sec  for  incident  waves  with  periods  from  1 to  3, 

3 to  7,  and.  7 to  10  min,  respectively.  Resonant  peaks  were  then  further 
defined  by  calculation  intervals  of  0.5  to  2 sec  and  sometimes  0.1  to 
0.25  sec.  If  a resonant  peak  appeared  for  one  harbor  configuration  and 
not  another  or  was  larger  for  one  configuration,  an  effort  was  usually 
made  to  locate  the  peak  and  resolve  it.  Occasionally,  it  was  difficult 
to  locate  a very  sharp  peak  or  it  was  not  considered  crucial  to  resolve 
peaks  for  some  gages.  Ship  motion  problems,  for  example,  are  of  more 
concern  in  the  East  Channel  than  at  gage  locations  LA-1  and  LA-5.  Very 
shatq  resonant  peaks  also  are  probably  not  important  in  the  prototype 
since  only  a small  amount  of  energy  in  the  incoming  wave  spectrum  would 
be  amplified. 

50.  Plates  1+3-89  show  frequency  response  plots  of  wave-height 
amplification  for  10  gages  located  along  the  landfill  areas  of  the  mod- 
ified harbor.  Plots  for  gages  LA-6  through  LA-10  display  peaks  at  90, 
157,  and  258  sec.  LA-11  has  an  additional  peak  at  116  sec,  and  LA-12 
has  a broad  peak  at  186  sec  and  no  noticeable  peak  at  258  sec.  Plots 
for  gages  LA-13,  LA-lU,  and  LA-15  display  no  significant  resonant  peaks. 

51.  Ship  motion  problems  due  to  very  long  wave  excitation  are 
due,  of  course,  to  horizontal  water  movement  and  not  to  vertical  water 
movement.  Waves  with  periods  from  1 to  10  min  have  amplitudes  of  only 
several  tenths  of  a foot  at  most;  therefore,  the  heave  of  a ship  induced 
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by  waves  is  insignificant.  Horizontal  motions  of  sway  and  surge  may  be 
large,  however.  Plates  90-170  are  frequency  response  plots  of  normal- 
ized maximum  current  velocity  (NMCV).  The  NMCV  is  a maximum  horizontal 
velocity  over  one  wave  period,  and  it  is  normalized  so  that  multiplica- 
tion by  the  incident  wave  amplitude  in  feet  yields  velocity  in  units  of 
feet  per  second.  Incident  waves  in  the  prototype  may  have  amplitude 
orders  of  magnitude  less  than  1 ft,  and  the  water  velocities  they  pro- 
duce are  correspondingly  smaller. 

52.  Plates  90-99  show  the  NMCV  as  a function  of  wave  period  for 
prototype  gages  1 and  2.  Velocities  are  small  at  gage  1,  of  course, 
since  it  is  at  the  end  of  the  East  Channel  where  antinodes  usually  de- 
velop when  the  channel  is  excited.  Plots  of  the  NMCV  at  gage  2 show 
several  resonant  peaks  for  both  existing  and  modified  conditions. 

Gage  2 was  located  at  an  oil  terminal,  and  reports  of  heavy  ship  motion 
were  received  from  tankers  several  times  during  the  prototype  data  col- 
lection of  Reference  10. 

53.  The  NMCV  plotted  as  a function  of  wave  period  is  shown  in 
Plates  100-10b  for  the  sta  2 LNG  slip;  this  gage  location  is  near  the 
mouth  of  the  LNG  slip  as  shown  in  Plate  2.  A NMCV  plot  is  not  shown 
for  the  sta  1 LNG  slip  because  the  velocity  plot  is  very  similar  to  the 
plot  for  the  sta  2 LNG  slip,  but  the  velocities  are  lower.  Plates  100, 
102,  and  103  show  resonant  peaks  at  110,  216,  and  258  sec.  Spatial 
velocity  plots  of  the  NMCV  are  shown  for  the  110-  and  258-sec  plot  in 
Plates  17b  and  188,  respectively. 

5^.  Plates  105-123  show  the  NMCV  as  a function  of  wave  period 
for  gage  locations  sta  LA-1,  LA-3,  LA-b,  and  LA-5.  LA-1  was  located  at 
the  end  of  the  West  Channel  in  Los  Angeles  Harbor  and  consequently 
there  are  no  large  peaks  in  the  NMCV  plots  shown  in  Plates  105-109. 
Gages  LA-3  and  LA-b  were  both  in  the  East  Channel.  Plates  110-119 
show  plots  of  the  NMCV  for  these  gages.  Velocities  are  quite  large  for 
the  b-  to  8-min  period  range.  The  two  major  resonant  peaks  from  b to 
6 min  are  larger-  for  the  modified  harbor  conditions  than  for  the  exist- 
ing conditii  ns;  however,  these  peaks  are  fairly  sharp.  The  response  of 
the  East  Channel  to  a spectrum  of  incident  wave  periods  may  be  larger 
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for  existing  conditions  as  a result  of  the  larger  current  velocities 
for  wave  periods  in  a range  between  the  two  peaks.  Velocities  are 
larger  for  existing  conditions  for  the  6-  to  8-min  period  range. 

Plates  120-123  indicate  that  the  NMCV  values  at  gage  LA-5  are  quite 
similar  for  existing  and  modified  conditions. 

55-  Gages  LA-6  through  LA-11  are  within  the  proposed  general 
cargo  terminal  to  be  created  by  the  Terminal  Island  landfill  shown  in 
Plate  2.  The  NMCV  plots  shown  in  Plates  121-153  for  these  gages  all 
display  resonant  peaks  at  90,  157,  2l6,  238,  258,  and  370  sec  and  a 
broad  peak  from  305  to  350  sec.  The  magnitudes  of  these  peaks  change 
from  gage  to  gage.  • Whether  or  not  these  resonant  peaks  would  induce 
ship  motion  problems  depends  upon  response  characteristics  of  the  type 
of  ship  moored  at  these  locations  in  the  terminal  area. 

56.  No  large  resonant  peaks  are  seen  in  the  NMCV  plots  for  gages 
LA-12  through  LA-15  shown  in  Plates  I5U-I7O.  Plate  2 shows  the  loca- 
tions of  these  gages. 

57.  Plates  171-198  are  plots  of  contours  of  wave-height  amplifi- 
cation and  plots  of  normalized  maximum  current  velocity  for  many  of  the 
resonant  peaks  shown  in  Plates  11-170.  The  velocities  are  represented 
by  arrows  whose  centers  lie  at  element  centroids.  Water  particles  move 
horizontally  back  and  forth  in  a direction  parallel  to  the  arrows.  The 
amplitude  of  the  incident  waves  is  assumed  to  equal  1 ft.  The  veloc- 
ities are  in  units  of  feet  per  second  for  this  amplitude  incident  wave. 
Incident  waves  in  the  prototype  may  have  amplitude  orders  of  magnitude 
less  than  1 ft.  Plates  171-198  basically  illustrate  patterns  of  wave- 
height  amplification  and  current  velocity,  that  is,  they  show  how  the 
harbor  area  is  oscillating. 
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PART  V:  CONCLUSIONS 

58.  The  good  agreement  between  prototype  measurements  in  the  Los 
Angeles  and  Long  Beach  Harbors  complex  and  calculations  of  the  hybrid 
finite  element  numerical  model  confirms  the  reliability  of  this  numeri- 
cal model  as  a viable  engineering  tool  for  studying  resonant  oscilla- 
tions of  complex  harbors.  Amplitudes  of  resonant  peaks  predicted  by 
the  numerical  model  will  always  be  too  large  since  the  model  does  not 
include  any  dissipative  effects  other  than  radiation  dissipation.  How- 
ever, the  model  can  predict  relative  severities  (except  that  friction 
is  likely  to  increase  with  increasing  amplitude)  and  is,  therefore,  ex- 
tremely useful  in  comparing  the  oscillation  characteristics  for  exist- 
ing conditions  with  those  for  proposed  modifications.  The  rapid  calcu- 
lation time  of  the  numerical  model  makes  it  feasible  to  apply  the  model 
to  large  harbors  and  broad  incident  wave  period  bands.  This  numerical 
model  can  be  used  as  an  aid  to  hydraulic  model  tests  since  it  pinpoints 
the  periods  of  resonant  peaks  and  shows  oscillation  patterns  which  are 
useful  in  optimizing  wave  gage  placement  to  measure  extreme  conditions. 

59-  The  numerical  model  calculations  indicate  that  resonant  os- 
cillations in  existing  harbor  facilities  in  the  Los  Angeles  and  Long 
Beach  Harbors  complex  will  not  be  substantially  altered  by  the  proposed 
modifications  of  the  harbors  complex.  Resonant  oscillations  in  Long 
Beach  Harbor  are,  in  general,  somewhat  smaller  for  the  n.  dit'ied  con- 
ditions of  the  harbor  complex  than  for  existing  conditions  (Refer- 
ence 11).  Several  peaks  are  larger  for  the  modified  conditi  ns;  how- 
ever, they  are  not  significantly  larger.  Some  resonant  peaks  in  exist- 
ing harbor  facilities  of  Los  Angeles  Harbor  are  increased  and  others 
decreased  by  the  proposed  modifications  of  the  harbors  complex.  How- 
ever, most  of  the  peaks  which  are  increased  are  sharp  peaks.  Sharj 
resonant  peaks  probably  are  not  important  in  the  prototype  since  only 
a small  amount  of  energy  in  the  incoming  wave  spectrum  would  be  con- 
tained in  the  narrow  period  range  of  the  sharp  peak.  The  broad  peaks 
displayed  in  the  frequency  response  plots  presented  in  this  report  for 
gages  in  existing  facilities  are  usually  similar  for  existing  and 
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modified  harbor  conditions.  The  greatest  difference  in  broad  peaks  of 
the  wave  amplification  factor  and  the  KMCV  for  existing  and  modified 
harbor  conditions  occurs  over  the  period  range  from  I4OO-5OO  sec  in  the 
East  Channel.  The  broad  peaks  for  the  modified  harbor  conditions  are 
substantially  less,  than  those  for  existing  conditions  over  this  period 
range  in  the  East  Channel. 

60.  The  largest  resonant  peaks  of  the  NMCV  in  the  LNG  slip 
occurred  at  incident  wave  periods  of  110,  157,  and  258  sec.  The  largest 
pieaks  of  the  KMCV  in  the  general  cargo  terminal  created  by  the  Terminal 
Island  landfill  occurred  at  wave  periods  of  90,  157,  and  258  sec.  Spa- 
tial plots,  of  the  NMCV  and  contour  plots  of  the  wave-height  amplifica- 
tion  factor  for  the  incident  wave  periods  of  90,  110,  157,  and  258  sec 
are  shown  in  the  text  in  order  to  completely  define  the  resonant  modes. 
This  information  would  be  required  as  input  to  a numerical  model  which 

rulated  shi]  res]  nse.  seaward  side  of  the  Terminal  Island  land- 

fill lid  not  exhibit  significant  resonant  oscillations  for  incident 
waves  with  periods  from  i to  lo  min. 

61.  Detailed  informal  i.  t.  n moored  ship  response  as  a function 

of  incident  wave  am;  l i • ai>  and  ;■  rl  d,  the  incident  wave  spectrum,  and 
the  resonant  response  f the  hurl*  .r.  ••  implex  to  incident  wave  excita- 
tion must  be  known  bet1  re  iefiniti  ■'  >nclus  L >ns  can  be  made  on  the  pre- 
cise degrees  of  shi]  mot!  n Li  Angt  les  and  Long  Beach  Harbors.  Nu- 

merical model  calculation.:  pc  . -nte  1 in  this  report  indicate  that  the 
resoi  . responses  of  existing  hard  r facilities  in  the  Los  Angeles  and 
Long  Beach  Harbors  complex  w u.  1 pr  t ably  riot,  be  changed  enough  to  cause 
an  increase  in  the  severity  : resonant  oscillations.  Whether  or  not 
ships  using  facilities  it.  ' r. ■ pr  p,  • i harbor  expansion  of  the  Los 
Angeles  Harbor  w il  1 exp<  rienc.  r...  riant  iscillation  difficulties  de- 
pend. upon  the  e harao t • ri s'  i the  ships  themselves.  Problems  would 

deve|o|  sly  if  the  r at  . . ft . ci  l 1 at.  ions  of  the  ship  were  close 

in  peri  d f r ••  •ilia*:  r.s  (redid'  d by  the  numerical  model  for 

basin  areas.. 


* 
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Table  1 


Es 

timates  of  Wave 

Periods  and  Spectral 

Energy 

Density  of  Peak  Energies  for  Gage 

_1 

Starting  Time  of 

Ranges  of  Wave  Periods 

Time  History 

1+80-360 

170-95 

80-60 

35-25 

71-166-1800 

387.02* 

133-93 

69.52 

30.21 

739- 7** 

231.6 

88.09 

. 

71-168-0300 

387.02 

133-93 

71.55 

33.51 

386. 1+ 

171+.1+ 

61+.  98 

9 . 666 

71-289-1600 

387.02 

120.77 

67.61 

29.85 

279-9 

180.1+ 

78.82 

10.70 

71-30^-1200 

387.02 

li+1.65 

78.39 

29 . 81+ 

1+51-5 

39-53 

13.1+1+ 

8.717 

71-31U-1500 

387.02 

120.77 

69.52 

29.RU 

533-9 

182.9 

86.19 

10.52 

71-316-0900 

387.02 

115.11 

67.61 

30..  ! 

1+08.5 

82.1+5 

1+6.07 

6. 1+63 

71-321,-0000 

387.02 

133.93 

69.52 

29.81+ 

1+83.1+ 

68.19 

28.66 

5.21+1+ 

71-329-0300 

387.02 

133.93 

73.69 

29.81+ 

660.9 

282.5 

1+8.70 

13.27 

71-331-1500 

387.02 

115.11 

65.80 

29.11+ 

1+28.1 

39.72 

lit.  81 

8.086 

72-155-0600 

387.02 

133.93 

67.61 

29.1+8 

311.1+ 

97.05 

66.26 

11.01 

* Wave  periods  in  seconds.  0 0 

**  Spectral  energy  densities  in  10  1 ft‘  -sec. 
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APPENDIX  A:  NOTATION 


a Boundary  of  region 

a^  Amplitude  of  incident  wave 

a Amplitude  of  transmitted  wave 

A Region  inside  harbor 

b Incident  wave  amplitude 

0 

b(w)  Amplitude  of  frequency  component  w 

d Average  rock  diameter 

2 

g Acceleration  due  to  gravity,  32.2  ft /sec 
G Element  slope  matrix 

h Water  depth,  ft 

H Hankel  function  of  the  first  kind  of  order 
n 

1 Imaginary  number 

k Wave  number 

n Integer 

n Unit  normal  vector  outward  from  region  A 
a 

N Interpolation  function 

r Spherical  coordinate,  ft 

R Real  number 

e 

t Time , sec 

T Wave  period 

u Velocity  in  x-direction,  ft/sec 

U Total  horizontal  velocity,  ft/sec 

v Velocity  in  y-direction,  ft/sec 

w Angular  frequency,  radians/sec 

x Cartesian  coordinate,  ft 

y Cartesian  coordinate,  ft 

a Unknown  coefficient 

n 

Y Coefficient 

A Area  of  element 

V Gradient  operator,  ft-"*" 

0 Spherical  coordinate,  radians 

A Wavelength  of  incident  waves 
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Response  of  harbor 
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Total  velocity  potential,  ft  /sec 

2 

Total  velocity  potential  evaluated  on  boundary  a , ft  /sec 

2 

Velocity  potential  of  incident  wave,  ft  /sec 

2 

Far  field  velocity  potential,  ft  /sec 

2 

Scattered  wave  velocity  potential,  ft  /sec 
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